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Abstract. We present an analysis of all BeppoSAX ob- 
servations of NGC 4151. This source was observed 5 times 
from 1996 to 2001 with durations ranging from a day to 
four days. The intrinsic continuum (described as a cut-off 
power law) , is absorbed at low energies by a complex sys- 
tem: a cold patchy absorber plus a warm uniform screen 
photoionized by the central continuum. We find that this 
"dual absorber" is the main driver of the observed vari- 
ability, up to a factor of eight, at 3 keV. In particular 
the covering fraction of the cold absorber changes on time 
scales of the order of a day, supporting its association with 
the Broad Line Region. The column density of the warm 
gas varies on a longer time scale (months to year). Some 
of the small amplitude spectral variability above 10 keV 
can be explained with an intrinsic variation (with varia- 
tion of the photon index AF ~ 0.2). The flux below 1 keV 
remains constant confirming an extended origin. Its spec- 
trum is reproduced by a combination of a thermal compo- 
nent (with temperature kT = 0.15 keV) and a power law 
with the same slope as the intrinsic continuum but with 
an intensity a few per cent. A Compton reflection com- 
ponent is significantly detected in 1996 (averaged value of 
ri/27r ~ 0.4, with fl the solid angle covered by the refiect- 
ing medium), with intensity decreasing on time scale of 
year, and it desappears in 2000 and 2001. The long time 
scale of variations argues for an association with an opti- 
cally thick torus at a distance of few light years. An iron 
line was detected in all spectra. Its energy is consistent 
with fluorescence by cold iron. We find that the line is 
variable. Its behaviour is reproduced by a variable com- 
ponent proportional to the level of the reflection flux plus 
a constant component. The flux of the latter is consistent 
with the extended line emission observed by Chandra. We 
conclude that the first component is likely arising from the 
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torus and the second is produced in the extended Narrow 
Line Region. 



Key virords: Galaxies: Seyfert; X-rays: galaxies; Galaxies: 
individual: NGC 4151 



1. Introduction 

The Seyfert galaxy NGC 4151 is one of the brightest 
AGN. It has been considered for years the prototype of 
its class and, as such, has been extensively studied at all 
wavelengths (see Ulrich 2000 for a review). In X-rays its 
spectrum is the most complex observed so far in AGN, 
being characterized by narrow and broad spectral fea- 
tures from soft to hard X-rays (e.g. IPerola et al. 19861 
IZdziarski et al. 19961 IZdziarski et al. 2002p . The central 
X-ray continuum is absorbed below a few keV by a 
(complex) absorber with column density TVh ~ 10^^ 
cm~^ associated with the Broad Line Region (BLR, 
[Reichert, Mushotzky fc Holt 1986| . This absorpt ion m 
NGC 4151 allows to reveal the presence of soft X-ray 
components below 1 keV, associated with a scattering 
component and a thermal plasma with a low tempera- 
ture of ^ 0.1 keV (|PeFola et al. 1986MPounds et al. 19861 
IWeaver et al. 1994b|) as those observed in type 2 objects 
(e.g. lAntonucci 19931 IMatt et al. 1997p . 

At higher energies NGC 4151 shows an iron line 
which remains remarkably constant notwithstanding 
large variations of the continuum (jPerola et al. 198*61 
ISchurch et al. 20^3|) . Observations by Chan- 
dra ( [Ogle et al. 2000D have indeed shown that most 
of the line is produced in the Narrow Line Region 
(NLR). The iron line site of Seyfert 1 galaxies has 
often been associated in the past with an accretion 
disk (e.g. ITanaka et al. 1995j IFabian et al. 2000p . XMM- 
Newton and Chandra observations have shown that a 
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narrow line component is a common feature in Sy Is 
(Reeves et al. 2001; Pounds et al. 2001; Matt et al. 2001; 
Kaspi et al. 2001). Further insight on the nature of 
the reprocessing medium should also become available 
through measurements of the reflection component. 
Recent INTEGRAL observation (|Beckmann et al. 2005]) 
described the hard X-ray/soft-gamma spectrum of 
NGC 4151 with a Compton continuum from hot electrons 
{kTe ~100 keV) in an optically thick (r=1.3) corona plus 
a reflected component from a cold material subtending a 
solid angle of il/2TT ~ 0.7. 

NGC 4151 has been one of the first Seyfert 1 galaxy 
with a well established evidence of changes of the spectral 
index correlated with the luminosity IjPerola et al. 19861 



Yaqoob et al. 1993). With the increased bandwidth and 



sensitivity of present instrumentation, this feature appears 
to be present in other Seyfert 1 galaxies (e.g. IH 0419-577, 
Page et al. 2002; NGC 3783, De Rosa et al. 2002; MCG- 
6-30-15, Vaughan & Edelson 2001; NGC 5548, Petrucci 
et al. 2000, Nicastro et al. 2000; NGC 7469, Nandra 
et al. 2000; IC 4329A, Done et al. 2000; Mkn 509, De 
Rosa et al. 2004, IZdziarski et al. 2003^ . In summary, m 
the broad X-ray range from 0.1 to 200 keV, NGC 4151 
represents a case study of the common properties (absorp- 
tion/scattering medium) of type 1 and type 2 objects, of 
AGN environment from light-day to kpc, and of investi- 
gation of the properties of the intrinsic continuum. 

In this paper we present the results of the complete 
set of observations by BeppoSAX performed in 1996, 1999, 
2000 and 2001. The observations and data reduction are 
described in Section [2l A model-independent variability 
study is presented in Section [3l The broad-band spectral 
analysis is reported in Section |4] and discussed in Section 
[5l Our conclusions are drawn in Section [6l 

2. BeppoSAX observations and data reduction 

NGC 4151 was observed by Bep- 
poSAX dPiro, Scarsi &: Butler 19951 IBoella et al. 1997|) 
several times: in July 1996, December 1996, January 
1999, December 2000 and December 2001. The log of the 
observations is given in Tabled] 

The first observation, part of the Science Verification 
Phase (SVP), was originally planned to last for about 3 
days. However the pointing was interrupted for 1 day and 
then restarted 2 days after. In the observations in 1996 
the three MEGS units, the HPGSPC and the PDS were 
operated. The LEGS was always on, with the exception 
of the second part of the observation in July 1996. From 
January 1999 onward two MEGS units, the LEGS and the 
PDS were operating. 

The data of all instruments were reduced using stan- 



dard procedures (Fiore, Guainazzi & Grandi 1999). A BL 



Lac object is located ~ 5' away from NGC 4151. The 
count rates are ~ 2 per cent and ^ 20 per cent of that of 
NGC 4151 for the MEGS and LEGS respectively. For the 



MEGS we adopted the standard extraction circle of 4', in 
which the contribution of the BL Lac is negligible ( ^ 0.4 
per cent). In the extraction of spectra and light curves of 
LEGS we excluded a circle of 2' of radius centered on the 
BL Lac object, while adopting the standard source extrac- 
tion radius of 8'. We estimate a residual contribution of 
the BL Lac of « 5 — 10 per cent, comparable to the statis- 
tical errors. The PDS spectra filtered either with fixed or 
variable rise time thresholds did not show any significant 
difference. The latter method is more suited to sources 
fainter than NGC 4151, so we have adopted data derived 
with the first method. 

In Figure [T] we show the light curves of the LEGS (0.1- 
2 keV), MEGS (2-10 keV) and PDS (13-100 keV). The 
MEGS (2-10 keV) counts rate during the July 1996 ob- 
servation changed by a factor of two on time scales of 
one day, from a low state (hereafter J96L) to a high state 
(J96II). The intensity during the December 1996 obser- 
vation was constant (D96), and quite similar to the J96L 
state, then low. Also the observation in December 2000 
(DOO) was constant. In January 1999 the average flux was 
at an intermediate level between the Low and High states 
defined above. During the observation the source flux ex- 
hibited a slow increase in the first part of the observation 
(J99L1), reached a peak (J99H) and then started to de- 
crease (J99L2). The time scale was ^ 0.5 days, with an 
amplitude of - 1.8 (max to min). The PDS (13-100 keV) 
light curve follows the same pattern of variability of the 
MEGS, but with a reduced amplitude: for example the 
ratio high to low states in July 1996 is 2.00 ± 0.01 for 
the MEGS and 1.22 ± 0.03 for the PDS (Table [J). On 
the reverse the LEGS (0.1-2 keV) count rate, with the 
exception of December 2001 (DOl), does not show signif- 
icant variability: the small difference (< 20 per cent, see 
Table [T]) of level between different observations is con- 
sistent with the effect of grid obscuration of the source 
(Fiore, Guainazzi & Grandi 1999), placed in different po- 
sitions in the three observations. The different behaviour 
in December 2001 will be discussed in Section [46l Consid- 
ering the substantial amount of spectral variability shown 
by this object, we have produced separate spectra for each 
of the above mentioned states, in which (2-10 keV) vari- 
ations are ^ 25 per cent. 

Hereafter errors and upper limits on spectral parame- 
ters correspond to Ax^ = 2.7, i.e. 90 per cent confi- 
dence level for a single parameter of interest. Spectral 
power indexes refer to photon distribution. In spectral 
fitting the normahzations of LEGS, HPGSPC and PDS 
relative to the MEGS have been left free to vary within 
the current ranges of uncertainty: LEGS/MEGS: 0.7-1.0; 
HPGSPC/MECS=0.95-1.1; PDS/MECS=0.77-0.95. 

3. Spectral ratios 

The spectral ratios derived comparing the spectral states 
observed during our campaign are presented in Figure 121 



Table 1. Log of the observations 



Date 


+1 

''exp 


LEGS (0.1-2 keV) 
cnt rate^ 


''exp 


MEGS (2-10 keV) 
cnt rate^ 


HPGSPC (8-20 keV) 

tl^p cnt rate^ 


PDS (13-200 keV) 
texp cnt rate® 


6-10/7/96 (J96L) 


6 


0.065 ± 0.003 


56 


1.145 ± 0.005^ 


26 


1.86 ± 0.06 


26 


2.36 ± 0.02 


(J96H) 






15 


2.29 ±0.01^ 


7 


2.7 ±0.1 


8 


2.90 ± 0.04 


4-7/12/96 (D96) 


13 


0.051 ± 0.003 


48 


1.102 ±0.005^ 


23 


1.9 ±0.06 


19 


2.28 ± 0.03 


4-6/1/99 (J99L1) 


14 


0.053 ± 0.002 


36 


0.937 ± 0.005'' 






16 


2.55 ± 0.04 


(J99H) 


10 


0.057 ± 0.003 


27 


1.238 ± 0.007'* 






13 


2.95 ± 0.04 


(J99L2) 


10 


0.059 ± 0.002 


20 


1.015 ± 0.007'' 






9 


2.60 ± 0.04 


22-23/12/2000 (DOG) 


4 


0.035 ± 0.003 


19 


0.470 ± 0.005" 






9 


1.01 ±0.04 


18-21/12/2001 (DDI) 


43 


0.145 ±0.002 


114 


1.929 ± 0.004* 






53 


2.90 ± 0.01 



Note: ^ net exp. time in ks; ^ For the MEGS we indicate the total rate of all MEGS unit operating; ^ 3 MEGS units; * 2 MEGS units; ® s"^ per half 
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Fig. 1. Light curve of NGC4151 in cts/s. We note that from 1999 onward (lower panels) two MEGS units were 
operating, while in 1996 (upper panels) they were three. 



This allows us to derive a model-independent description 
on the spectral variability of the source. We note the fol- 
lowing: 

1. Below 1 keV the flux remains constant, notwithstand- 
ing the large variations showed at E ^ 2 keV. The low 
energy spectrum is thus to be dominated by constant 
component (s) (jPerola et al. 19861 IPounds et al. 19861 
IZdziarski et al. 1996P . This means that some absorber 
has to suppress any variable intrinsic component at 
low energies. Around 0.6-0.7 keV the ratio DOl/DOO 
shows evidence of a changing of the opacities of the 
warm gas. We will discuss this effect in Section l4?2l 

2. The 2-10 keV range is characterized by the largest 
variability, up to a factor of ~ 8 in the 2-3 keV range. 
The spectrum gets steeper when the flux increases 
(J96H/J96L, J99H/J99L1, DOl/DOO, J99H/J96L). 
This happens both on short (days) and long time scales 
(years). This behaviour has been attributed to intrinsic 
spectral variations correlated with the luminosity on 



time scale as short as few hours (Yaqoob et al. 1993 
IPerola et al. 1986]) . and to variations of the absorber 
on month-year time scales. Warwick et al. (1995) and 



Weaver et al. (1994b) instead attributed all variations 
to changes of absorber structure. We will assess this is- 
sue through a broad band spectral fitting (see Section 

3. Above 10 keV two kind of behaviour are apparent, 
a) The spectral ratios relative to observations per- 
formed within « day are all consistent with a constant 
(J96H/J96L: xl = 1-2, J99L2/J99L1: xl = 0.6 and 
J99H/J99L1: xl = 1-2). b) The spectral ratios relative 
to observations performed over time scales of months- 
years are not (J96H/J99H: xl = 2.0 , and DOl/DOO: 
xl = 1.8, and J99H/J96L: xl = 2.6). This is a new 
result that, as it will be shown in Section 14.31 bears 
important implications on the origin of the reflection 
component. 



In the next sections we will present a detailed quanti- 
tative analysis of these results, through detailed spectral 
deconvolution. 
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Fig. 2. Spectral ratio between different states. 



4. Broad-band (0.1 200 keV) spectral modelling 

On the basis of previous observations and models de- 
velopped in previous published work, to fit the Bep- 
poSAX broad-band spectra (with all the instruments in 
the energy range specified in Table [T]), we adopted a spec- 
tral template model (baseline morfe/BLM), which includes 
the following components: 

A. The intrinsic continuum is described by a power-law 
with an e- folding energy Ec- 

B. A Compton reflection component (PEXRAV model 
in XSPEC, Magdziarz &: Zdziarski 1995| , with the cosine 
of inclination angle of the reflector fixed to 0.95. We 
assume elemental abundances from Anders & Grevesse 
(1989). 

C. A narrow Ka iron line modeled with a gaussian 
profile with the intrinsic width set to 10 eV. 

D. The soft X-ray spectrum {E < 2 keV) is a combi- 
nation of a thermal component model (MEKAL) with the 
temperature fixed to 0.15 keV and a scattering component 
(i.e. a power-law with the same slope of the intrinsic hard 
continuum but with a different normalization). 

E. The complex absorber is reproduced with a dual 
absorber model. A fraction fcov of the source is covered by 
a cold column density A^hcov An additional uniform and 
photoionized gas is responsible for additional absorption 
of the spectrum at low energies and for a midly ionized 
Fe edge at 7.3-7.5 keV (Fe VII-FeVIII) (ABSORI model in 
XSPEC IDone et al. 1992)) . 

This photoionized gas is characterized by a tempera- 
ture T—3 X 10* K, a column density A^Hwarm and a ion- 
ization parameter ^ = Lion/nR^, where Lion is the source 
luminosity between 5 eV and 300 keV, n is the hydrogen 
number density of the gas and R its distance from the cen- 



tral ionizing source. The Fe abundances in the cold and 
warm absorbers and in the reflection component are tied 
together and left free during the fit. 

Our model assumes that only the e-folded power-law 
with reflection (i.e. the intrinsic continuum), is subject 
to complex absorption. Additional absorption through the 
Galaxy is applied to all the emission components (A^Hgai = 
2 X 10^0 cm-2). 

The BLM model is shown in Figure [3l In the upper 
panel in Figure the spectrum with model of J96L is plot- 
ted with all the different spectral components while in the 
lower panel two spectral states with a flux variation of a 
factor of ~ 3 (see Table [2]) are shown: December 2000 and 
December 2001. We applied this template to all the Bep- 
poSAX spectra in the different flux level states (see Table 
[1] and light curves in Figure [T|) . The best fit parameters 
and the reduced values are shown in Tabled while the 
data/model ratio for each spectrum is plotted in Figure 
m In the case of July 1996 high state (when LEGS data 
are not available) the parameters of the soft components 
(thermal emission and scattering component) are frozen 
to those derived for the spectrum in July 1996 low state. 

The broad-band baseline model described above pro- 
vides a good fit to all the spectra of NGC 4151 (Tabled]). 
In the following sections we discuss the various spectral 
components and their variability. 

4-.1. The origin of the spectral variability in the 2-10 keV 
range and the intrinsic power law 

The source has always shown - in all the observations per- 
formed so far and including ours - a systematic trend: the 
spectrum gets softer when the flux increases. The origin 
of this behaviour is one of the key issue in this source. 



Table 2. Broad band (1.8-200 keV) fits with our BLM. Eq is the e-folding energy of the intrinsic continuum, Aic the normahzzation of the intrinsic 
continuum at 1 keV, Area is the strength of the Compton reflection hump obtained as fl/27r x ^ic, Ifb is the flux of the iron hne at 6.4 keV, A^hcov and 
-^Hm'arm axc thc cohimn densities of the cold partially covered and warm gas respectively, /cov is the covered fraction of the cold gas, ^ is the ionization 
parameter of the warm gas, E'Mthermai is the emission measure of the thermal gas - difined by / neUHclV, where rie is the electron density, and nu is the 
hydrogen density - assuming the distance to the source of 13 Mpc, and Agcatt is the 1 keV flux of soft scattered component. 



Period F ^E^ '^Ajc ^^rcfl '^^Fc (■^/^q)fc "^-^Hcov /cov '^-^Hwarm ^ £, ^_EMtherm ^^scatt ^-fT-fokcV ^ fi)!?- 200keV Xi 
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Note: (1^ e-folding energy in keV; in 10"^ keV"! cm-2 s'l; in 10"'' cm-2 s"! at the line; in 10^^ cm-2; (s) in erg cm s'l; in IQ-i" erg cni-2 s'l; ^'^^ in 
10®^ cm"^ ; in 10"^ keV"i cm s ; * These parameters are frozen to their values of best fit computing the errors. 
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Fig. 3. BeppoSAX broad-band spectrum EF{E). In the upper panel we plot the spectrum in July 
1996 low state with all the spectral components of the BLM described in the text: in green (solid 
line) the total spectrum, in red (dot line) the scattered component, in grey (dash-dot line) the 
thermal soft comnonent. in blue and maeenta fdash-dot-dot lines^ the transmitted comnonents 
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Fig. 4. The Data/model ratios for the basehne model. 



The broad band fit (presented in Table [21) shows that which can be related to the 2-10 keV spectral variability, 
the only spectral parameters showing significant changes, 
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are those of the absorbing media. Indeed, a decrement of 
the effective absorption (i.e. of A^h and/or /cov) results in 
a softer spectrum, as appears when comparing the spectral 
ratios presented in Section [3] with the best fit parameters 
in Table[2] For example, the softening observed from J96L 
to J96H (Figure [2|) is reproduced by decreasing the col- 
umn densities of the two absorbers and the covering frac- 
tion (Tabled]). On the contrary the intrinsic photon index 
does not show significant evidence of variation around the 
average value of F « 1.6. A variation AF w 0.2 would 
still be consistent with the errors, but it could not itself 
account for the large spectral changes observed in this 
range. While most of the spectral variability observed in 
the 2-10 keV range is due to variation of the absorbers, 
a minor effect of intrinsic variations can not be excluded. 
This is demostrate by fitting the spectra with the intrinsic 
photon index given by the F2_iokeV vs F relation found 
by Perola et al. (1986): F = 2.18 + 0.012 x -F2-iokoV with 
^2-iokeV in 10~^^ erg cm~^ s~^. The results of these fits 
are reported in Table[3]and show that in all spectra a good 
xi is obtained. 



4-. 2. The complex absorber 

With the exception of the 2001 observation, the cold ab- 
sorber is characterized by A^hcov = (1 — 3) x 10^"^ cm~^ 
and /cov = 0.3 — 0.7, similar to that found by Weaver 
et al. (1994b), who adopted the same model of this anal- 
ysis on ASCA data. The model employed for the warm 
gas with ^ 1-16 erg cms^^ indicates an iron edge 
at energy around 7.3-7.6 keV, consistently with other 
measurements ( Yaqoob et al. 1993 IWeaver et al. 1994a| 
IWeaver et al. 1994bl ISchurch et al. 2003j) . With the ex- 
ception of DOl, the iron abundance is w 1.5 times greater 
than the solar value. We will show later how DOl is pecu- 
liar in respect to several points. 

With regard to the variability of absorbers, we note 
the following. Due to the fit with such a complex model, 
the errors on the absorbers parameters are wide. However 
we note that the column density of the warm absorber 
appears to show the most significant variations in respect 
to those of the cold absorber. We have therefore chosen to 
fix the column density of the cold absorber to the average 
value TVhcov = 15 x 10^^ cm~^. This is consistent, within 
the errors, with all the best fit values except DOl. In ad- 
dition we have fixed the spectral index, as from the above 
discussion, to the value 1.6. The result and all the best fit 
parameters are reported in Table [H The is still accept- 
able. The column density of the warm absorber changes 
over a time scale of month-years, while it is consistent 
with being constant over a day time scale. Specifically, in 
July 1996, A^Hwarm « 3.5 x 1022 cm-2, in D96 to DOO 
Nnvjarm ~ 7 X 1022 cm~2. Qu the contrary the covering 
fraction of the cold absorber shows variability in a day 
time scale. For example changing from 0.73 to 0.48 in the 
low and high state of July 1996. 



The DOl spectrum appears to be rather peculiar with 
respect to the other states. We already noted that in 
this state the iron abundance is about 5 times that so- 
lar. In DOl the softening is reproduced with a substan- 
tial decrease of the thickness of both absorbers, with 

AThcov = 3.5 X 1022 cm-2, iVHwarm = 0.9 X 1022 Cm-2. 

The covering fraction was also very low, /cov ~ 0.3. Ba- 
sically, in the state observed by BeppoSAX in December 
2001, most of the bare nuclear continuum came into view 
(Figure [31 lower panel). The DOl spectrum is peculiar 
also for the presence of an absorption feature at 9 keV. 
This state of the source is described and discussed in Piro 
et al. (2005). 



4-. 3. The reflection component 

The Compton reflection component is significantly de- 
tected in 1996 (both in July and December, with a mean 
value of the absolute normalization ^refl = (2.2 ± 0.3) x 
10^2 ph cm-2 g-i ]jgV-i)j it is only marginally detected 
in 1999 {AreR = (0.8 ± 0.3) x 10-2 ph cm-2 g-i keV-^), 
while it is consistent with zero in 2000 and 2001 (see Ta- 
ble |4|). This result remains true even when the photon 
index is not fixed (see Table [2] and Table [3]). Thus we find 
that the reflection component is variable on time scales 
of years, while there is not strong evidence of variability 
on the (intra) day time scales that characterize the direct 
continuum. We can actually exclude variations correlated 
with those of the intrinsic continuum on such short time 
scales. 



4-. 4- High- energy spectral variability 



In Section 14.21 we have concluded that the spectral vari- 
ability observed in the 2-10 keV range is dominated by 
variations of the absorber. At higher energies, although 
the effect of absorption is negligible, the assessment is also 
not straightforward due to the presence of the reflection 
component. For example, a constant absolute reflection 
component, added to a steep (constant slope) variable con- 
tinuum (jSchurch fc Warwick 2002)) . may also qualitatively 
mimic the F vs i^2-iokeV correlation. It is then particu- 
lary relevant to discuss the high energy spectra in 2000 
and 2001 when the reflection intensity is consistent with 
zero and then any intrinsic spectral variation can be ob- 
served without ambiguity. Between 10 keV and 200 keV 
the spectral ratio DOl/DOO is not consistent with a con- 
stant {xt — 1.8 per 15 degree of freedom), but shows a de- 
creasing trend that can be reproduced with a change of the 
intrinsic slope between the two states of about 0.17 ±0.06 
(with the assumption of a constant high energy cut-off). 
This result is that expected on the basis of the correlation 
F vs F2-iokcV, where F2-iokcV is changing by a factor of 
three between DOO and DOl. 



Table 3. As Table [5] with the intrinsic photon index F frozen following the F2_iokoV vs F relation found bv lPerolaet al. 19861 



Period r 'E, 2^IC ^ A,,a '^h. (^/^0)fc ''A^Hcov ,fcov ''iVnwarm ^SMthorm 



J96L 


1.42 


68 ±5 


3.0; 


yvi:2 

-0.2 


1 1+"'* 


3.9; 


fO.5 
-0.5 


9 0-1-0.4 
^•■-'-0.3 


20.1+^'"7 


66+" "^ 

"•""-0.06 


5.ill;l 


10* 


8.11^;? 


9 9+0. « 
^•^-0.5 


2.04 


96.4/78 


J96H 


1.54 




5.2: 


hO.7 
-0.5 


1 6+"'^ 


2.9; 


fO.9 
-0.9 


2.9+lt 


10.4+?-^ 


o.5o«:S« 


2.6+?-? 


16* 






3.01 


73.2/68 


D96 


1.42 




2.8: 


1-0.3 
-0.2 


1 Q-fO.4 


4.1; 


fO.5 
-0.5 




18.lli°- 


0.54«:1° 


8.31?;^ 


8.3* 


10.0+^'? 


2.0l«;« 


2.01 


65.6/78 


J99L1 


1.51 


10011° 


4.6: 


1-0.4 
-0.2 


> 0.5 


4.2; 


fO.7 
-0.8 




19.6+11 


"•"■J-o.os 


7.91°:^ 


8.7* 


i3.oi?;i 


i-8lg:l 


2.72 


61.9/73 


J99H 


1.59 


122+1^ 


6.5: 


[-0.5 
-0.4 


Q n + 0.6 


2.9; 


fO.9 
-0.9 


9 1-1-0.5 
^•-■--0.4 


17.3+t^ 


"•00-0.07 


6.3lJi 


2.7* 




1 r+O.S 


3.38 


86.2/73 


J99L2 


1.47 


loollt 


3.9: 


[-0.3 
-0.3 


^•^-0.2 


4.2; 


fO.9 
-0.9 




10.9l|;^ 


"•O'-0.28 


D.0_2.o 


10.9* 


5.11^5:? 


9 f,+0.7 
^•"-0.5 


2.44 


73.1/73 


DOO 


1.30 


llOl'as 


i.o: 


[-0.1 
-0.1 


< 0.1 


1.9: 


f 0.5 
-0.5 




2.111:^ 




1.6l?i 


10* 


9 q+4.4 
^•"J-1.6 


1 q+0.3 
^^•'3-0.3 


1.01 


74.8/73 


DOl 


1.50 


1221^ 


4.7: 


[-0.1 
-0.1 


< 0.2 


8+0-4 
2-^-0.4 


n 1 -f2.2 
'•J^-2.1 


q 1+0.9 

■"■■^-o.e 


41+0 09 
"•^1-0.07 


o.sig;^ 


0.3* 


2-8l?:^ 


9 q + 0.5 

^•'3-0.9 


2.66 


70.3/73 



Note: e-folding energy in keV; in 10"^ keV"^ cm"^ s'^ in 10~* cm"^ s"^ at the line; in 10^^ cm"^; in erg cm s"^ in 10"^° erg cm"^ s" 
10^^ cm~^; in 10~^ keV~^ cm~^ s~^;* These parameters are frozen to their values of best fit computing the errors. 
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Table 4. As Table [21 without less important component, with r=1.6 and the cold column density A^hcov = 15 x 
lO^^cm-^. 



Period ^E, MiC 'A.efl '/fc /cov ''iVHwarm xl 



J96L 




q q+U.iJ 

■^■•J-o.i 


q q+0.6 
"-'••-'-0.3 


4 2+"-* 

^•^-0.4 


n 7q+0.U3 
"• ' "^-0.04 




8* 


104/79 


J96H 


98tlt 


"•"-0.8 


1 9+0.6 
.^•^-0.6 


n 4+1.1 
■^•^-0.5 


0.48t°°^ 




16* 


73/69 


D96 


124+1^ 


q c;+0.2 
'-'■^-0.2 


9 9+0.5 
^•^-0.3 


^•-■--0.4 






9.9* 


76/79 


J99L1 


135tf5 




^^•"-0.4 


^•^-0.7 


n 70+0.04 
"• '■-'-0.06 


6.9t;« 


8.6* 


73/74 


J99H 


127+^4 




0^6^o'5 


2 g+0.9 
^•''-0.9 


"•""-0.05 


"•"-0.7 


3.2* 


87/74 


J99L2 


144+^° 


r 1+0.4 

^■J^-0.4 


0.8«-J 


4 9+0.9 
^•^-0.9 


n 70+0.07 
"• ''-'-0.16 


74+3.0 


9.7* 


74/74 


DOO 


/1 99+3965 
'±^^_225 


1 S+0.1 


< 0.5 


^•"-0.6 


"•'J*-0.28 


6.7+?-^ 


14* 


71/74 


DOl 


162tn 


"■"-0.2 


O+0.3 
"•'-'-0.2 


^•"-0.4 


r, 90+0.02 
"•^^-0.03 


1 0+0.2 

^^•^-0.1 


2.8* 


107/74 



Note: keV; in 10"^ koV"^ cm"^ s"^ in 10~* cm"^ s"^ at the line; in 10^^ cm"^ 



December 2001: absorbed power law model Decembre 2001: dual absorber model 




6 ° 6 

channel energy (keV) channel energy (keV) 



Fig. 5. In the left panel we show the iron line profile when the continuum is reproduced with a simple absorbed power 
law in 3-10 keV. Clear positive residuals are evident between 5-6 keV. When the continuum is correctly modeled with 
a dual absorber model (our BLM), the residuals disappear, on the right panel, and a narrow gaussian profile becomes 
a good description to the iron line. The broadening of the profile in right panel is due to the MECS resolution (8 per 
cent at 6 keV). 



4-. 5. The iron line 

We recall that, in our modelling, we have adopted a 
narrow gaussian profile for the line without any further 
component, following the results of Chandra and XMM- 
Newton ( |Ogle et al. ^OOOl ISchurch et al. 2003p . Because 
the issue of the existence and origin of broad components 
in Seyfert galaxies is still under debate, it is relevant to un- 
derline how this is dependent upon a good modelling of the 
continuum, available thanks to the BeppoSAX broad-band 
data. In Figure [H] we show the residuals when the contin- 
uum is modelled by a simple uniform absorbed power law 
and the line shows evidence of a "red wing" . However, 
when the continuum is properly described with a dual ab- 
sorber, the "red wing" - which is in fact due to the covered 
component - disappears. 

The flux of the iron line in our spectra shows evidence 
of variability {x^/dof = 19.6/7 for a constant, chance 
probability P=0.0065), but these variations are not cor- 
related with the intensity of the intrinsic continuum (cor- 



relation coefficient r=0.099). Interestingly the flux of the 
line is lower in DOO and DOl when the reflection compo- 
nent disappears. To explore this issue we plot in Figure 
|6] the line vs reflection intensity. A linear flt, plotted in 
the flgure, is acceptable {x^ /dof — 7.8/6, P=0.253) with 
/ = (2.5±0.3) X 10-^-h(0.57±0.17) x Ur'^Ar^a, I is in ph 
cm^^ s^^and Arefi in ph cm~^ s^^ keV^^. This flt shows 
that when the reflection component is null, the line flux 
does not go to zero, identifying the presence of a constant 
line component in addition to the one that correlates with 
the reflection. While it is not statistically signiflcant, it is 
nonetheless interesting to note that for large value of the 
reflection normalization, the line appears to saturate. This 
point will be discussed later. 

4. 6. The soft components 

In the soft X-ray band the spectrum is well described by 
the combination of a power law continuum and a thermal 
plasma with fcr=0.15 keV (see Table [2|). In Table [5] we 
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Fig. 6. The flux of the Iron line is plotted as a function 
of the Compton reflection component. We plot also the 
result of a linear fit 

Table 5. Flux in 0.1-1 keV energy range. We report the 
total flux and that of the thermal and power law compo- 
nent. All the values are corrected for Galactic absorbtion. 



Period 


jptot 


ptherm 


piscatt 


J96L 


6.6 


1.4 


5.1 


D96 


6.6 


1.2 


5.4 


J99L1 


6.7 


2.3 


4.4 


J99H 


6.8 


3.0 


3.8 


J99L2 


5.4 


0.6 


4.8 


DOO 


4.7 


0.7 


4.0 


DOl 


8.1 


1.6 


5.5 



Note: Flux corrected for Galactic absorption, in 10~ erg 
cm~^ . The J96H state in not reported because the 
missing of soft X-ray data from LEGS (see Section [2}. 



report the flux of these two components in the 0.1-1 keV 
range. The total flux in this range is is dominated (« 80 
per cent) by a constant power law with i^(0.1 — 1 keV) « 
5 X 10""'^^ erg cm^^ s^^, the remaining being attributed to 
the thermal component. We have computed the expected 
flux of the variable absorbed continuum in this energy 
range, flnding that it is at most 1.5 per cent of the total. 
The only exception is the soft spectrum DOl where there 
is a contribution of the intrinsic continuum in 0.1-1 keV 
due to the effect of low absorption (see Table ^ . 

5. Discussion 

5.1. The primary emission from the innermost region 

The intrinsic continuum of NGC 4151 is well described 
by a power law with an exponential cut-off. The val- 
ues of the high energy cut-off we found in the different 
flux level states are all consistent with the values ob- 
served in the Sy Is by BeppoSAX (jPerola et al. 2002^ . 



This, together with the measure of the Compton re- 
flection features that perfectly match with the averaged 
values observed in Sy Is (jNandra fc Pounds 1994^ . sug- 
gests that NGC 4151 is an intrinsically average Seyfert 
1 galaxy (jZdziarski et al. 2002[) . The phenomenological 
model we employed to fit the intrinsic continuum (i.e. 
power-law with high energy cut-off), is a good descrip- 
tion of a two-phase model involving a hot corona emit- 
ting medium-hard X-rays by Comptonization and a cold 
optically thick layer (the disk) that provides the soft 
photons to Comptonization (jHaardt fc Maraschi 1991] 
[Haardt, Maraschi fc Ghisellini 1997D . While the dual ab- 
sorber is the main driver of the observed variability in 
2-10 keV, we find small amplitude of intrinsic variations 
above 10 keV, that can be due to a variation of AF w 0.2, 
according to the the correlation F vs i^2-iokoV (Perola 
et al. 1986). Assuming a distance of 13 Mpc for the source, 
and a black hole mass of 4x 10^ Mq we find that the ratio 
-'^o"i^20okcv/-'^E between the unabsorbed luminosity and 
the Eddington luminosity is in the range 0.002-0.006 in 
the lowest (DOO) and highest (J99H) state respectively. 
The accretion is taking place in a sub Eddington regime. 

5.2. The narrow iron line and the reprocessing region: 
NLR and molecular torus 

The energy of the narrow line as determined from the aver- 
age of all BeppoSAX observations is (6.40 ±0.06) keV, i.e. 
consistent with fluorescence by cold iron. During our ob- 
servations, the line flux shows evidence of variations (see 
Table ll}. We find that the flux of the line is dependent 
upon the normalization of the reflection component, and 
that it attains a non-zero value (/ = (2.5 ± 0.3) x 10^'' ph 
cm~^ s~^) when the reflection is absent (in Figure [6] we 
plot the flux of the iron line as a function of the reflection 
amplitude) . This suggests that the line is made up by two 
components. Wc first discuss the origin of the component 
at zero reflection. Our result indicates that this compo- 
nent is consistent with being constant, thus suggesting an 
extended origin. In addition, it should be produced by 
an optically thin gas. The flux we derive for this compo- 
nent is remarkably similar to the spatially resolved line 
emission region observed by Chandra ( |Ogle et al. 2000| 
/ = (1.8 ±0.2) X 10"** ph cm~^ s~^). This region is coinci- 
dent with the NLR, also in agreement with our conclusion 
of an optically thin gas. 

The equivalent width of iron line in the case of an opti- 
cally thin medium illuminated by an isotropic continuum 
is given by (e.g llnoue 19891 IPiro 19931) : 

EW^im^^{^),.^eY (1) 

where Ai7/47r is the net solid angle covered by the 
medium (equal to the angular extension of the region 
multiplied by /cov)- We can estimate the column 
density of the NLR A^h = fRn, by adopting the typi- 
cal values derived from optical observations of the size 
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i? « 1 kpc, density (n ^ 200 - 1000 cm- ^) and filling fac- 
tor (/ « 10"''; e.g. lRobinson et al. 1994p . By substituting 
these values in eq.(l) we derive an EW ~ 0.1 — 0.3 eV 
that is 2 orders of magnitude lower than observed (100- 
200 eV see Table [2]). One possibility is that the ioniz- 
ing flux is emitted anisotropically. Pension et al. (1990) 
showed that the ionizing flux needed to produce the ob- 
served optical and UV line emission had to be about 10 
times stronger than that observed. Explanations for this 
excess include an anisotropic radiation field. But even in 
such a case, the expected EW would be lower of the ob- 
served value by a factor of 10. The presence of a system 
of clouds with higher density (n « 10^ cm~'^) can pro- 
vide a more appealing explanation. In this hypothesis the 
column density of NLR would be A^h 3 x 10^"^ cm^^, 
then from eq. 1 we get an equivalent width EW ^ 200- 
300 eV, that is good agreement with the values we found. 
Spectral diagnostics of narrow lines in a sample of Seyfert 
nuclei that includes NGC 4151, reveals the presence of low 
density {n ^ 10** cm""^) and high (n ^ 10^ cm~'^) density 
clouds (jStasihska 1984p . Moreover, imaging of the NLR of 
NGC 4151 with HST shows the presence of regions whose 
density is likely to be enhanced by the shock compression 
of the ejected radio material ( |Winge et al. 1997| . 

Let us now discuss the line component that is corre- 
lated with the reflection normalization, characterized by 
/ = (0.57 ± 0.17) X lO-^Arefl ph cm-2 s-i(see Figure 
[6]). The line produced by reprocessing by a cold thick 
medium has an intensity / « 8 x 10~^G'y [(Z/Z0)Fe]^rofl 
ph cm^^ s~^where G is a function depending on the ge- 
ometry, and is equal to « 1 (substantially independent 
on the inclination angle) for simple configurations like a 
planar disk ( |Matt, Perola fc Piro 199T| Y[{Z/Zq)f^] de- 
scribes the non-linear dependence of the line EW upon the 
iron abundance. For {Z/ZQ)-pe = 2, the case of NGC 4151, 
it is y = 1.6 (jPiro 1993|) . and the line flux should be 
I Ki 1.2 X 10~^GArefl ph cm~^ s~^, remarkably similar 
to the observed one for G « 0.5. Values of G less than 
unity are expected in a geometrically and optically thick 
torus, where the reprocessing takes place mostly in the 
inner wall, i.e. in a cylindrical or conical geometry. This 
is due to the "secondary illumination" effect, produced by 
photons emerging from the wall at angles that intercepts 
the wall again. This effect enhances the efficiency of repro- 
cessing and does it more for the photons of the reprocessed 
continuum than for those in the line, because a 6.4 keV 
photon has a high probability to be absorbed. For exam- 
ple, if the disk is in a conical configuration with an opening 
angle of 45°, G « 0.5 ( |Matt, Perola fc Piro 19911 ). 

We now discuss the origin of the variability of the re- 
flection component and the associated iron line. This could 
be due either to a variation of the illuminating continuum 
or to a change of the geometry of the reprocessing region. 
We recall that the intensity of the reflection decreases on 
a time scale of years by a factor of two or more from 1996 
to 2000-2001 (see Table H]). The size of the reprocessing 
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Fig. 7. Light curve in 1.5-10 keV of NGC 4151 from ASM 
aboard RXTE . Each bin is the flux averaged on 60 days. 

region is therefore to be of the order of light-years. The 
intensity of the reflection component thus traces the level 
of the illuminating continuum averaged and delayed over a 
time scale of « year. We have therefore looked into the his- 
torical light curve of NGC 4151, to search for evidence of 
long time scale variability. In Figure [7] we report the light 
curve as observed by the ASM aboard of RXTE. Each 
point is the average flux over 2 months. The source shows 
a systematic decrease in its intensity by almost a factor 
of two from 1996 to 2000, consistent with the observed 
reflection variability. 

The most obvious candidate for the reprocessing re- 
gion is a thick torus whose inner walls are at « light years 
away from the central source. Taking then into account 
the line and reprocessing intensities we estimat^ that the 
distance of the inner region of the torus should be « 3 — 4 
times its height. The ionizing continuum would therefore 
emerge from the torus within a cone of opening angle of 
« 70°, in agreement with the estimation based on the ex- 
tended Narrow Line Region made by (jPedlar et al. 19831 
IRobinson et al. 1994p . Finally we note that the variabil- 
ity of the reprocessing component could alternatively be 
explained by a change of the height of the torus. 

5.3. The X-ray absorbers and the BLR 

BeppoSAX observations of NGC 4151 show the presence 
of a complex absorber system. In particular the intrinsic 
continuum is absorbed by a cold gas that partially obscure 
the central source and by a second uniform warmer gas. 
Our analysis shows evidence of variations in the status of 

^ A precise determination of the geometry of the reproces- 
sor goes however beyond the scope of this paper. It would re- 
quire detailed computations, that should take into account the 
non-linear effect of secondary illumination, the solid angle sub- 
tended by the inner surface of the torus towards the X-ray 
source and their exposed area towards the observer. 
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the cold absorber on time scales of the order of the day 
and of the warm gas on longer time scale. Variations of the 
covered component can be produced if the cold medium 
is made up of clouds passing in the line of sight towards 



emission. We reproduced the soft X-ray spectrum with 
two components. One, dominated by emission lines, which 
is well fitted by a thermal plasma with kT « 0.15 keV 
with abundances comparable to solar. This is in agreement 



the central source (Reichert, Mushotzky & Holt 1986). If 
we identify these clouds with the BLR, from the typi- 
cal velocity of ~ 10^ km s~^, we derive that the size of 
the central source R « 10^'* cm. This would correspond, 
for i? w 10 Schwarzschild radii, to a black hole mass 
of 4 X 10^ Mq, consistent with independent estimations 
(|Clavel et al. JMH [Wamlel, Peterson fc Malkan 1999| ). If 
the cloud velocity is due to the orbital motion, we de- 
rive that the distance of the clouds is a few light days, 
consistent with the size of the inner BLR in NGC 4151 
(|Clavel et al. 19871) . 

In our fit we have associated the warm stage with 
a uniform screen covering the source. In doing so we 
are guided by the arguments put forward by Warwick 
et al. (1995) to solve the discrepancy between the X-ray 
and UV measurements of absorption. Upper limits on UV 
extinction (jKriss et al. 1992]) correspond to a column of 
< 8 X 10^° cm~^ assuming a standard gas-to-dust ratio. 
This compares to X-ray columns of 10^^ — 10^^ cm~^. In 
our case, the level of photoionization is such that a sub- 
stantial fraction of light elements as H, He and C are ion- 
ized fe.g. lKrolik fc Kallman 1984^ . and therefore their ab- 
sorption which is predominant in the UV range is reduced 
to a level compatible with UV measurements. We there- 
fore associate the uniform absorber to the warm medium, 
while the patchy absorber is in a cold (neutral) stage. 

Let us now characterize this warm medium we associ- 
ated to the UV line absorbing system (jKriss et al. 1992]) . 
From the relationship ^ — L/nR^ « 10 erg cm s""'^, with 
Lo.i-20okcV ~ 3 X 10*'^ erg s~^ and the lower limit on the 
density n ^ 10^-^ cm~^ derived by UV observations, we 
obtain that the size of this region should be i? ;S 10^^ cm. 
This is consistent with the upper limit derived from the 
X-ray variability {R ^ 10^^ cm) and the requirement that 
the medium be external to the BLR (i? ^ 10^^ cm). In 
conclusion, our data are fully consistent with the scenario 
in which the ionized X-ray absorber coincides with the 
system responsible of UV absorption lines. This system 
is composed by many clumps of dense material covering 
the broad line region, and produced by an outflow of the 
broad-line clouds (jKriss et al. 1992^ . 

The substantial decrease of the warm and cold absorb- 
ing gas in DOl together with the evidence, in the same 
spectrum, of an absorption feature at ^ 8.5-9 keV sug- 
gests the presence of a multi-phase ionized absorber. This 
topic is widely discussed in Piro et al. (2005). 



use components: scattering and 



5.4- Low energy 
thermal 



The low energy emission in NGC 4151 remains constant on 
time scales of the order of years, indicating an extended 



with Chandra observations ( Ogle et al. 2000 ) that demon- 
strated that most of the soft X-ray flux is extended and 
dominated by lines from a photoionized and coUisionally 
ionized plasma. The other component is a power law with 
a normalization « 3 — 10 per cent (see Table 3) of that of 
the hard power law, which is likely to be produced by elec- 
tron scattering of the central continuum, as in the case of 
Seyfert 2 galaxies (jAntonucci 1993|IMatt et al. 1997^ . The 
fraction of this component in respect to the continuum is 
strongly dependent from the intrinsic flux. We expect the 
scatterd fraction increase when the source flux decreases 
if the nuclear continuum is scattered into our line of sight. 
Weaver et al. (1994a) in ASCA observation of NGC 4151 
found a fraction scattered of 3-4 per cent for the high state 
and 5-6 per cent for the low state. In our BeppoSAX ob- 
servations the upper value of the ratio Agcatt/^ic ~ 10 per 
cent is obtained in the DOO spectrum when the continuum 
level was lower in respect to the other states. 

6. Conclusions 

In this paper we presented the broad-band (0.1-200 keV) 
BeppoSAX spectrum of NGC 4151. The source was ob- 
served 5 times from 1996 to 2001 with durations ranging 
from a day to four days. We find that the spectrum is com- 
plex characterized by several components. The continuum 
emission is well reproduced by a cut-off power law. 

At energies less than 5 keV the continuum is strongly 
obscured by a dual absorber: a cold component, that we 
associated with the BLR clouds orbiting around a cen- 
tral source with a black hole with a mass of 4 x lO^M©, 
characterized by iVn « 1.5 x lO^'^ cm^^, a covering frac- 
tion /cov=0.3-0.7 and variable on time scales of days to 
years, and a second warmer uniform screen (with iVu = 
(1 — 8) X 10^^ cm"^, variable on time scales of months- 
years), photoionized by the continuum emission and that 
we argue to be coincident with the UV absorbing region. 

The low energies spectrum is reproduced by a combi- 
nation of a two components, which are not absorbed by 
the cold and warm gas: a thermal one and a power law 
with the same slope of the intrinsic continuum but with 
the intensity which is a few per cent. 

During our observations the source shows strong time 
variability energy dependent in time scales ranging from 
hours to years. In 0.1-2 keV it did not show any variability, 
while in 2-10 keV the variation was up to a factor of eight 
at 3 keV and three at 10 keV. In the PDS energy range 
(13-200 keV) the larger amplitude of variation was of a 
factor of three. We found that in 2-10 keV most of the 
spectral variability can be attributed to a variation of the 
absorbers. Instrinsic spectral changes can explain some of 
the smaller amplitude spectral variability above 10 keV. 
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Wc detected a Compton reflection component at high 
energy and the iron hne at 6.4 keV. The intensity of the 
reflection component is significantly different from zero in 
1996 while it disappears in 1999 and later (2000-2001). 
The size of this region, implied by the variability argu- 
ment, as well as its optical thickness argue for an associ- 
ation with a geometrically and optically thick torus sim- 
ilar to that of Seyfert 2 galaxies. The study of the his- 
torical light curve of NGC 4151 suggests that the reflec- 
tion component traces the level of the illuminating con- 
tinuum averaged and delayed over a time scale of year. 
We found that the iron line flux is well described by two 
components: one is correlated with the reflection inten- 
sity and is associated with the torus, the other is constant 
I = (2.5 ± 0.3) X 10~^ ph cm~^ s~^(in very good agree- 
ment with Chandra measurement). This component has to 
be produced in a diffuse optically thin g confirmed 
by Chandra observations that have shown that about 70 
per cent of the line is extended over a size of w 1 kpc 
coincident with the extended NLR. 
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